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ABSTRACT 

An analysis of IGR J 16207 -5 129 is presented based on observations taken with Suzaku. The data set 
represents ~80ks of effective exposure time in a broad energy range between 0.5 and 60keV, including 
unprecedented spectral sensitivity above 15keV. The average source spectrum is well described by an 
absorbed power law in which we measured a large intrinsic absorption of A^h = (16.2;^j j) x 10^~ cm"^. 
This confirms that IGR J 16207-5 129 belongs to the class of absorbed HMXBs. We were able to constrain 
the cutoff energy at 19^4keV which argues in favor of a neutron star as the primary. Our observation 
includes an epoch in which the source count rate is compatible with no flux suggesting a possible eclipse. 
We discuss the nature of this source in light of these and of other recent results. 

Subject headings: accretion, accretion disks ; gamma-rays: general ; stars: neutron ; supergiants ; X-rays: binaries ; 
X-rays: individual (IGR J 16207-5 129) 



Introduction 



Seven years after its launch, INTEGRAL dWinkler et al 



caUed "IGRs" (for INTEGRAL Gamma-Ray source^ 
now number over 300. While INTEGRA L'S gamma- 



2003h continues to discover sof t gamma-ray sour ces in 
the 20-lOOkeV energy range jBird et al.ll2010l) . So- 



ray imager ISGRI jLebrun et alj l2003h discovered 



these objects, it provides an error radius that is typi- 



for an updated list of 
|http://irfu.cea.fr/Sap/IGR-Sources/| 



IGRs, 



please 



Table 1 
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observation ID 


dates 


instrument 


exposure time 


name 




(MJD) 




(ks) 




402065010 


54499.823-54500.784 


HXD 


49.596 


Ol 


402065020 


54526.866-54527.750 


HXD 


28.331 


02 



XIS 32.678 



cally around 2^', i.e. too large to allow for the identi- 
fication of an optical/lR counterpart which would help 
classify the source into one of the many families of 
high-energy emitters. Therefore, follow-up observa- 
tions wdth^ocu^ingX^ray telescopes such as Chandra 
(e.g. Tomsick et al.ll2006t) . Swi ft (e.g. Rodriguez et al 



2010l) . and XMM-Newton (e.g. iBodaghee et ai1l2006l) 
are necessary to refine the X-ray position and to extend 
the spectrum below ISGRJ's ~20keV lower limit. 

One such source, IGR J 16207 -5 129, was discov- 
ered in the first mosaic images constructed from 



among the growing class of absorbed SGXBs. Other 
spectral features typical of this class were unveiled 
such as a Ka-emission line at 6.39 + 0.03 keV from 
iron fluorescence, and soft excess emission at ~ 2 keV 
consistent with the re processing of X-rays in the wind 
jTomsick et al.ll2009l) . 

Suzaku is playing an increasing role in provid- 
ing spectral and timing follow-up observations of 



IGRs (e.g IBozzo et al.i l2008t iMorris et al.1 l2009t 



Paizis et al.ll2009l) . Recently, we asked Suzaku to tar- 



INTEGRAL-lSGm scans of the Galactic Plane jWalter et al 
2004t iTomsick et al.l 120041) . A follow-up observa- 
tion with Chandra allowed Tomsick et al.l ( 20061) to 
refine the source position to within an error radius 
of 0.6" which excludes all but a single optical/IR 
counterpart (2MASS J 16204627 -5 130060 = USNO- 
Bl.O 0384-0560875) whose spectral energy distribu- 
tion is that o f a hot and massive late O or early B- 
type star (e.g.lTomsick et al. 2006HMasetti et al.jl2006l: 



Negueruela & Schurchl2007tlRahoui et alJ2008l:lNespoh et al 
20081) . The svstem is therefore classified as a super- 



giant HMXB (SGXB). 

The X-ray characteristics of IGR J 16207 -5 129 
lend further credence to this classification. Although 
its emission shows variability by at least an order 
of magnitude, IGR J16207-5129 is persistent having 
been detected in short and long-expo sure images taken 



get IGR J16207-5129. Here, we present the results 
from the analysis of these data including: an X-ray 
spectrum in the 0.5-lOkeV range that confirms the 
obscured environment of the X-ray source; a spectrum 
with unprecedented sensitivity between 15 and 60keV 
allowing us to constrain the cutoff energy; and a light 
curve that includes a possible eclipse. The key steps 
of the analysis are presented in Section[2] followed by 
the pertinent results from timing and spectral studies 
in Section[3] Our conclusions and the implications of 
our findings are discussed in Section|4] 

2. Observations 

Launched in 2005, the Suzaku space telescope fea- 
tures 2 co-aligned instruments working in tandem to 
observe the hard X-ray sky (0 .2-600 keV): the X- 
ray Imaging Spectrometer (XIS: iKov ama et ak 20071) 
and t he Hard X-ray Detector (HXD: Takahas hi et al 



20061 2007, 20101; Krivonos et al. 


[2007 


) and other X- 


ray telescopes (e.g. Tomsick et al. 


I2OO6 


; Masetti et al. 



An absorbed power 
law fits the 0. 3-lOkeV Chandra sp ectrum well with 
F = 0.5:^°-^ jTomsick et al.1 l2006b . In the INTE- 
GRAL 20-50-keV band, the spectral slope steepens 
to F - 1.9 + 0.5 possibly indicating a high-energy 
cutoff jTomsick et alj I2OO6I) . Using XMM-Newton, 
Tomsick et al. ( 20091) measured a column density of 
A^H = (ll-9!o5) ^ 10^^ cm^-. This large intrinsic 
column density firmly establishes IGRJ16207-5129 



20070 . XIS consists of 2 front-illuminated CCD cam- 
eras (XISO and XIS3), and one that is back-illuminated 
(XISl). HXD contains 2 detectors, PIN and GSO, 
which are mainly sensitive to energies below 60keV 
and above 40 ke V, respectively. 

Focused on IGR J16207-5129, observation ID 
402065010 (PI: J. Tomsick) took place on Febru- 
ary 3^, 2008. A technical glitch affected the on- 
board processor for XIS resulting in no data from 
the low-energy instrument. Our observing time was 
rescheduled for March 1-2, 2008, as observation ID 
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Fig. 1. — Image of IGR J 16207 -5 129 as captured with the XISl detector on Suzaku. Pixels have been grouped in 
blocks of 4 (bin: group4 in ds 9). The circular source and squared background selection regions are also shown. 




Fig. 2. — Significance image (from 3 to lOcr) 
in the 17-80 keV band of IGR J 16207 -5 129 from 
INTEGRAL-lSGRl. The image is in Galactic coor- 
dinates and the white box represents the FOV from 
HXD-PIN. 

402065020. This time, both instruments recorded data 
successfully. All observations were performed in the 
XIS -nominal pointing mode while using a quarter win- 
dow option in order to gain a higher time resolution. 
Table [Upresents the observation log. 

Data were reduced with v.6.7 of the HEAsoft soft- 
ware package following the procedures described in 
the Suzaku ABC Guide v.2. The task xispi con- 
verted the unfiltered XIS event files to pulse-invariant 



channels, which were then input into xisrepro to 
return cleaned event files. After removing bad events, 
we extracted the source spectrum from each XIS de- 
tector using a circular region of 120"-radius centered 
on the position given by Chandra (see Fig. [Til. A 
pair of 240"-wide squares positioned on a source- 
free region of the detector served as the background. 
Response matrices were generated with xisrmf gen 
and xissimar f gen taking into account these source 
and background extraction regions. To improve pho- 
ton statistics, we merged the source spectral files for 
the front-illuminated CCDs (XISO and X1S3), and we 
did the same for the background. A minimum of 150 
counts was allocated to each bin in the XIS spectra. 

We used the tuned (v.2.0) background file and 
the common GTl during spectral extraction of HXD- 
PIN. We then performed dead-time corrections with 
hxddtcor. We modeled the cosmic X-ray back- 
ground using the flat-field response file from January 
29, 2008, and then subtracted its contribution from the 
HXD-PIN background spectrum. We tested whether 
the HXD-PIN spectra from Ol and 02 were consistent 
in terms of power-law slope (F ~ 2) and normaliza- 
tion before they were summed into a single spectrum. 
Channels in this spectrum were paired up between 15 
and 20 keV, grouped by fours between 20 and 40 keV, 
and collected in a single bin between 40 and 60keV. 
Spectral bins above this were excluded. 

Source and background light curves were generated 
in 2 energy ranges, 0.5-lOkeV (XIS) and 15-60 keV 
(HXD-PIN), for an effective exposure time of ~80ks. 
The time resolution was set to 2 s for XIS, and 32 and 
64 s for HXD-PIN. 

Figure|2] shows that the millisecond pulsar named 
PSR J1617-5055 is near the edge of the nominal 
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Fig. 3.— Background-subtracted light curve of IGR J 16207 -5 129 from XIS (red: 0.5-lOkeV) and HXD (blue; 15- 
60keV). The upper panel displays observation Ol while the middle and lower panels show 02 (see Table[T]for start 
times). Each bin collects 160 s worth of data. The hardness ratio is defined as where H and S are the count 
rates in 0.5-lOkeV and 15-60 keV, respectively. The dashed line represents MJD 54527.382 which corresponds to the 
onset of a period of low activity. 



HXD-PIN field of view (FOV). In the ISGRI energy 
band in which this image was captured (17-80keV), 
the ms-pulsar has an average flux around a quarter 
that of IGRJ16207-5129. According to the Suzaku 
technical description, the HXD-PIN count rate from 
a Crab-like source would decrease by a factor of 20 



when offset at 30' from the center of the FOV. This 
implies that the contribution expected from this source 
to the IGR J 16207 -5 129 count rate is 1-2% on aver- 
age, and up to 10% if the flux from IGRJ16207-5129 
diminishes by an order of magnitude. 
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ity ex pected of persistent SGXBs ( IWalter & Zurita Heras 
20071) . In con trast, supergiant fast X- r ay transients 



Fig. 4.— Full spectrum of IGR J 16207 -5 129 cor- 
rected for the background and fit with an absorbed 
power law with an exponential cutoff whose param- 
eters are listed in Table|2] The data represent photon 
counts from XISl (red), XISO combined with XIS3 
(blue), and HXD-PIN (black). 
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Fig. 5. — Combined XISO and XIS3 spectrum of 
IGR J 16207 -5 129 in the 5-8 keV energy range. Each 
bin collects 150 counts. An absorbed power law was 
fit to the data (solid curve). An iron line at 6.4 keV is 
suggested by the residuals to the model displayed in 
the bottom panel. 

3. Results 

3.1. Timing Analysis 

FigureO presents the background-subtracted light 
curves in 160-s bins from both soft and hard in- 
struments. The figure shows that the emission from 
IGRJ16207-5129 varies within 1-2 orders of magni- 
tude on ks-timescales, i.e. within the range of variabil- 



(SFX Ts: e.g. iNegueruela eTaP l2006t ISguera et al. 



20061) spend most of their time in a quiescent or in 
termediate phase (Lx ~ lO^^-lO^^* erg s"') punctuated 
by flaring episodes lasting a few hours th at represent 
variations of 2-5 orders of magnitude (iSidoli et al 



20081). The average (unabsorbed) flux during the ob- 
servation is ~ 3 X 10"" ergcm"^ s"' in the 0.5-lOkeV 
range which is similar to the averag e flux found with 



range wnicn is similar to tne averag e tiux rouna witi 
XMM-Newton jTomsick et alj 20091). Adopting a dis 



tance of 6kpc (iNespoli et al 



20081) . the luminosity 



recorded in XIS 1 at the highest point is 5 x 10 erg s" 
(between 24 and 26 ks after the start of 02), while it 
is 2 X 10-''*ergs"' at its lowest point (between 62 and 
72 ks after flie start of 02). 

Intriguingly, X-rays from the source, and the vari- 
ability associated with that emission, are significantly 
diminished during the last ~30ks of the observation. 
The flux, notably in the soft band, is consistently be- 
low the average registered for the full 02 data set be- 
ginning at around 45 ks into 02, which corresponds to 
MJD 54527.382. The light curves for 02 were then 
accordingly split into the epochs 02a and 02b before 
and after this date. 

For XIS, the median flux in 02a is 1.32 cps (std. 
dev.: cr = 0.7). This is higher by at least a factor of 
10 than the median flux in 02b: 0.11 cps (cr = 0.15). 
This difference between 02a and 02b is also reflected 
in higher orders of the distributions of light-curve bins: 
02a was less peaked and more symmetric in its distri- 
bution (kurtosis k - -1.2, skewness s - -0.1) than 
was 02b {k = 2.2, s = 1.3). The HXD Hght curves 
also show evidence of a prolonged attenuation (by a 
factor of 5) of photon count rates during 02b, with 
clear statistical differences in the higher orders of the 
distribution. The distribution of light curve bins from 
HXD-PIN displayed a sharper peak with more symme- 
try during 02a {k = -0.6, s = 0.03) than during 02b 
(;t = 9.8, i = 2.1). 

We searched the light curves from XIS and HXD 
for periodicities on the order of a few tens of seconds to 
a few thousand seconds. A moderate signal appears at 
~970s in the soft and hard bands (significance > 5cr). 
While this is in the range of reported spin frequen- 
cies for X-ray pulsars accreting from the winds of su- 
pergi ant companions dLiu et al. 1 120061: iBodaghee et al. 



20071 and references therein), this also coincides with 
a sixth of the orbital period of the Suzaku satellite. The 
970-s signal is also present in the background hght 
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Table 2 

Spectral parameters from models fit to XIS 0-3 and HXD-PIN spectra of IGR J 16207-5 129 



Xspec model 


C 


N„ 
10^2 cm-2 


r 


Em 
keV 


keV 


£ (K») 
keV 


norm. 

10--' 


xl/dof 




01+02 


c*phabs*pow 

c*phabs*highecut*pow 

c*phabs*highecut*(pow+gauss) 


1 05+°" 
96+°-" 


16.2!?;« 


1 36+00' 

'■■'"-0.07 

1 09+0" 

'■"^-0.09 

92+"--'' 


6 4+"-'' 

"■^-0.5 
4 5+2.6 
^■-"-0.5 


19!^ 


6 42+0-03 


4.4 
2.7 
2.1 


1.38/204 
1.15/202 
1.05/199 


02 


c*phabs*pow 


A 7^+0.08 
"■ ' -^-0.08 


18.7!|;» 


I ^9+0.08 
-0.08 








4.1 


1.30/212 


02a 


c*phabs*pow 


75+010 


19 l+'o 


1 •!4+0.09 
'■■^^-0.09 








6.2 


0.91/182 


02b 


c*phabs*pow 


'—'-0.56 


19.3!t;? 


'■^°-0.38 








1.4 


0.51/64 



Note. — Eitoi-.s are quoted at 90% confidence. C repre.sents an instrumental cros.s-calibration coefficient which is fixed to 1 for XIS 
and left free to vary for HXD-PIN. 



curve from HXD-PIN, as are fractional harmonics cor- 
responding to a fourth, and a third of the satellite's or- 
bital period. We conclude that this periodicity is due to 
aliasing of the satellite's orbit. It is worth noting that 
no periodic signal was detected i n a continuous ~50-k s 
observation with XMM-Newton dTomsick et al.ll2009l) . 



(iNespoh et al.ll2008h . the unabsorbed source luminos- 



3.2. Spectral Analysis 

Figure|4] presents the 0.5-60keV spectrum of 
IGRJ16207-5129 from all observations (Ol and 02 
combined). The source spectrum was initially fit with 
a power law tha t includes a photoe lectric absorption 
component with IWilms et al. I (l2000l) abundances, and 
an instrumental cross-calibration constant which was 
fixed at 1 for XIS, and left free for HXD-PIN. This 
fit gave a reduced value of 1.38 for 204 degrees of 
freedom (dof). Significant residuals remained above 
30 ke V, so we added an exponential cutoff which im- 
proved the quality of the fit (xl = 1.15 for 202 dof). 
The optimal fit (xl = 1.05 for 199 dof) was ob- 
tained by complementing this cutoff power-law model 
with a Gaussian line to account for FeKa emission 
at 6.4 keV. We expect to find this line in our spec- 
trum since it was d etected in a previou s observation 
with XMM-Newton dTomsick et al.ll2009b . The residu- 
als from the model without an iron line are shown in 
Fig.E] 

Parameters for these models are listed in Table|2] 
but we cite here specifically the column density of A^h 



r = 0.92 



+0.24 



and the cut- 



= (16.2+p X 10^^ c) 

off energy at 19^5 keV. Assuming a distance of 6kpc 



ity is 3 X 10 ergs in the 0.5-lOkeV range, which 
is twice the average luminosi ty measured with XMM- 
Newton dTomsick et al.ll2009l) . 



Visual inspection of the unbinned HXD-PIN spec- 
trum shows residuals to the power-law fit at ~22 keV, 
although there are no harmonically-spaced lines (specif- 
ically, at twice the energy of the purported fundamen- 
tal). We considered the possibility of a cyclotron line 
by adding a Gaussian absorption component (gabs in 
Xspec). This did not improve the quality of the fit 
compared to a power law. Furthermore, we could not 
simultaneously constrain the related optical depth, the 
line width, and line energy. The null hypothesis can 
not be excluded given an F-test probability approach- 
ing 80%. Maintaining the width of line at 3keV (a 
typical value found in the literature for cyclotron lines 
in other sources), we fixed the line energy to integer 
values ranging from 20 to 30keV, which we used to 
determine an upper limit of <3.4 on the optical depth. 

We sought to compare the spectrum of the source 
before and after MJD 54527.382 which corresponds to 
the onset of a noticeable change in the emission be- 
havior. Therefore, we split the XIS 0-3 and HXD- 
PIN spectra according to epochs 02a and 02b. Al- 
though the HXD-PIN recorded few photon counts dur- 
ing 02b, a careful grouping of the data enabled us to 
collect 4 usable bins between 15 and 60 keV. We then 
fit both spectra with absorbed power laws. The power- 
law slope and column density were consistent in both 



epochs: F = 1.34+|J|;; andA^n = (19!|)x 10 



22, 



for 



6 



02a; and T = 1.48!°-^ and A^h = (19!^) X lO^^cm-^ 
for 02b. The only significant difference between the 
spectra from 02a and 02b was in the normalization. 

Besides 02b, there were other instances in which 
the source count rate in XIS was low or close to 
zero, such as e.g. around 12ks and 36ks into 02a 
(see Fig.|3]l. Therefore, we separated the XIS spec- 
tra into low and high states corresponding to source 
count rates (prior to background subtraction) below 
and above 0.5 cps, respectively. The column den- 
sity from a power-law fit was slightly larger during 



4. Discussion 



the low state (A^h = (22.2+[ °) x lO^^cm^O, but re- 
mained statistically compatible with the high state (A^h 
= (19.4^12) X 10^^ cm"^). Once again, normalization 
of the spectral model remained the only significant dif- 
ference between the low and high states. 

Defining the hardness ratio as where S is the 
count rate in 0.5-lOkeV, and H is the count rate in 
15-60keV, we find that the source spectrum was soft 
at -1 during most of 02a, punctuated by 2 positive 
spikes during the brief episodes of low intensity at 12 
and 36 ks into 02. During 02b, the ratio tended to 
harder values, but since the dispersion is so large, we 
can not confidently state that the increase is signifi- 
cant. We redefine the hardness ratio, where S is now 
the count rate in 0.5-5 keV, and H is the count rate be- 
tween 5 and 10 keV. The boundary at 5 keV effectively 
splits the light curve into an equivalent number of to- 
tal counts, and ensures that variations in the column 
density will mainly affect the S -band. In this case, the 
hardness ratio during 02a varied between 0.5 and 1.0 
with hard spikes during the 2 brief low-flux episodes, 
and tended to be harder in 02b, but, once again, the 
uncertainties are large. 

According to the hardness ratios, the high-energy 
X-ray emission appears to be less affected than the 
low-energy emission during the low-intensity states of 
02a and during 02b. The nearby millisecond pul- 
sar PSRJ1617-5055 contributes 10% of the HXD- 
PIN emission during 02b which is not enough to ac- 
count for the observed change in the hardness ratio. 
We point out that the potential hardening observed in 
IGR J16207-5129 during 02b is similar to the ingress 
phase of the eclipsing SGXB 4U 1538-52 in which 
a non-negligible hard X-ray flux continues to be de- 
tected for 20 ks after the low -energy bandpass has 
reached a minimum count rate dRobba et al ] l200lh . 



4.1. The compact object and its environment 

The photoelectric absorption that we measured for 



IGRJ16207-5129 is A^h = (16.2+^^ x lO-^-^cm^O, 
which is 12-6 0% larger than the v alue found with 
XMM-Newton jTomsick et al.l l2009l) . and 3-7 times 
larger than the value reported with Chandra jTomsick et al 
20061) . suggesting a variable column density. The 
expected line-of-sight absorption (A^^) is among 
the highest in the Galaxy at (1.6-1.7jxl0^ ^ cm"^ 
dDickev & Lockmanll990HKalberla et al.li2005b . Sub- 
sequently, we infer that in addition to the interstel- 
lar medium that absorbs some of the radiation, there 
is matter obscuring the X-ray source that is local to 
the system. This situation is similar to, although not 
as extreme as, the ca se of IGR J 163 18 -4848 (e.g. 
FiUiatre & Chatvll2004 . 

We can estimate the density of material around 
SGXBs such as IGR J16207-5129 by integrating 
the wind de nsity along the lin e of sight to t he X- 



ray source ([Castor et al 
Prat et al.ll2008h : 



1975; Levine et al. 2004 



Pw 



Vwir) = Voo^l - 



(1) 



(2) 



Here, /?* is the radius of the stellar companion, Vco 
is the terminal velocity of the stellar wind, M* is the 
stellar mass-loss rate, and /3 is the wind velocity gradi- 
ent (a;0.8-1.2, we will assume /3 - 1). We adopt stan- 
dard values for B 1 la supergiant stars as provided by 
Crowther et al.l (120061) . specifically Voo = 725 km s ', 
and M,^ = lO^^Mgy"'. The radius of the companion 
star is set to R^, - 20 Rq. The compact object (CO) is 
considered to be a neutron star at the canonical mass 
of 1.4 Mq. Circular orbits are assumed with r - 2/?*. 
The binary is viewed edge-on, and the CO is at inferior 
conjunction. We also introduce the mass of the hydro- 
gen atom (niu) and the mean molecular weight of the 
accreted matter (/y ~ 0.6 for Solar material). 

Integrating Eqn.lT] we obtain: 



A^H 



A^ 



H+ 

Jr=2R, 

p dr 



r 

= A?S + 



(3) 
(4) 



7 



ln(2) 



(5) 



Inserting our assumed values into Eqn.[5] we derive 
a column density of 5 x 1Q~ cm"^. In order to be con- 
sistent with our observed value of A'h = 10^^ cm"^, 
we would need - 2 x lO^^Mgy"' which is a 
slightly larger, but still reasonable, stellar mass-loss 
rate. While this model can underestimate the true value 
of A^H for a given set of wind parameters, the derived 
column density should be higher if the system were 
viewed at superior conjunction (see Section l4~2t . A 
few caveats are in order. First, the intrinsic absorption 
in these sources can vary by an order of magnitude so 
a static model can not be expected to predict the full 
range of observed column densities. Second, the wind 
parameters that we assumed might not reflect the ac- 
tual conditions of this system. Third, we relied on an 
orbital geometry that is simplistic (e.g. the inclination 
angle is zero) since we know little about the binary's 
orbit. 

Iron fluorescence is another tell-tale sign of ma- 
terial near the high-energy emitter An iron line 
was detected for thi s source in XMM-Newton data 
(iTomsick et al. I l2009h . and our Suzaku spectrum in- 
cludes a possible iron line at 6.42+0.03 keV with an 
upper limit on the equivalent width at <9 1 eV, both 
of which are consistent with the corresponding values 
obtained with the XMM-Newton spectrum {E (Kff) = 
6.39+0.03 keV, EW =42+12eV). The line is produced 
by continuum X-rays fluorescing Fe tha t is restricted 
to neutral (Fe I) or ionized up to Fe XV (iHouse Il969t 
NagasJ 19891 and references therein). The line equiva- 
lent width combined with the column density are con- 
sistent with a model in which the fluorescence occurs 
within a spherically-symmetric shell of relatively cool 
matter (|Na gase"l989t lMat3l2002l:lKallman et alj|2004 : 
Drake et al.j i2008i) . Therefore, these measurements 
provide valuable clues to the properties of the matter 
around the CO. 

The nature of the CO in IGR J 16207 -5 129 remains 
unknown. Periodic modulations attributable to a neu- 
tron star's spin have not yet been detected from this 
source de spite long observations with XMM-Newton 
(ITomsick et al. 2009) and now with Suzaku. One pos- 
sible explanation for the non-detection of a pulsation 
is that the rotation and magnetic axes of the neutron 
star are nearly parallel. Another possibility is that 
an unfavorable geometric alignment causes the accre- 
tion beam to point in our direction throughout the 



rotation. As previously suggested by ITomsick et al. 
(l2009l) . IGR J 16207 -5 129 might have a very slowly- 
rotating neutron star, e.g. ~6ks s uch as in the mag ne- 
tar candidate IGR J16358-4726 dPatel et alj|2007l) . In 
our data, a signal at this frequency is too close to the 
orbital period of Suzaku to be unambiguous. 

Instead of a neutron star, the CO might be a radio- 
quiet black hole candidate. The hard X-ray emission 
from a black hole candidate generally has a cutoff 
energy situated above 60keV. In IGRJ16207-5129, 
the cutoff is at 20 ke V which suggests a neutron star 
primary (Nagase 1989 ). Soft excess emission is of- 



ten se en from accreting X-ray pulsars (iHickox et al 
2004 ), and this feature was detected in the XMM- 



Newto n spectrum of IGR J16207-5129 (ITomsick et al 
2009b . There were residuals around 2keV from the 



optimal fitting model to the XIS spectrum, but the ad- 
dition of a blackbody component did not improve the 
quahty of the fit in a significant way. 

Keep in mind that coherent pulsations have yet to be 
confirmed for many SGXBs including 
IGR J19140+0951, which shares many of the same 
properties as IGR J 16207 -5 129 such as the compan- 
ion's spectral type, a high i ntrinsic column density , 
and soft excess emission (e.g.lHannikainen et al .l2007t 
Prat et al.ll2008t iTorreion et al.ll2009l) . The lack of a 
measurable pulsation period or cyclotron absorption 
lines is common among SGXBs, and in particular 
among SFXTs, even though their spectral cutoff en- 
ergies argue in favor of a NS primary. 

4.2. Is IGRJ16207-5129 an eclipsing X-ray bi- 
nary? 

Our Suzaku observation ends with an epoch of 
~ 30 ks in which the source has an average count 
rate in the soft and hard bands that is 5-10 times 
lower than before this epoch, and up to ~20 times 
smaller at certain times. Whil e SGXBs typicall y 
vary in luminosity by a factor 20 dWalter et al.l(2006l) . 
the long duration of the inactivity in this source is 
unusual and therefore intriguing. During a 12-ks 
stretch of 02b, the unabsorbed 0.5-10-keV flux is 
4.3 X 10"'^ erg cm"^ s"' . This is lower than the lowest 
sustained flux (ks timescales) measured with XMM- 



Newton: 7.0x10'" ergcm'^s"' dTomsick et al1l2009l) 



In addition, the epoch of flu x minimum recorde d 
by XMM-Newton lasted ~3ks dTomsick et al.ll2009l) . 
which is 4 times shorter than the low-flux state that we 
observed with Suzaku. Note that the maximum flux 
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that we measured is consistent with the flux measured 
during the high-intensity state from the XMM-Newton 
observation: ~10"'"ergcm"^ s"'. 

There are several mechanisms that could suppress 
the source count rate from IGR J 1 6207-5 1 29 for a sus- 
tained period: an "off" state of the X-ray emitter; oc- 
culting clumps in the wind; eclipse of the shock region 
at the surface of the neutron star; an eclipse of the pri- 
mary itself; and hydrodynamic effects. In the next few 
paragraphs, we will discuss each of these possibilities 
and the extent to which they are supported by the data. 

The X-ray luminosity is driven by the accretion rate 
which is affected by, among other things, wind inho- 
mogeneities such as clumpy material, by photoioniza- 
tion, and by the magnetic field of the CO. In general, 
persistent emission with stochastic variability is ex- 
pected from such objects, but we point out that some 
SGXB s (e.g. VelaX-1), can enter short-lived inactive 
states (llnoue et al.lll984l) . The X-ray flux in this sys- 
tem decreases across all energies so varying absorp- 
tion columns are not responsible. The "off" states are 
thought to be due to the passage of the NS in a cav- 
ity of low-density stellar wind combined with an onset 
of the propeller regime which te mporarily inhibits ac 



cretio n onto the compact object dKrevkenbohm et al 



20081) . Typical timescales for these "off" states are on 
the order of a few hundred seconds, or at least an order 
of magnitude shorter than the period of inactivity that 
we observed in IGR J16207-5129. 

The X-ray source might be temporarily shielded 
from us by optically-thick blobs of the supergiant 
star's wind passing through our line of sight. We 
would anticipate an increase in A^h on short timescales 
concurrent with a decrease in the X-ray flux, poten- 
tially followed by a flaring episode depending on 
the proximity of the blob to the CO, i.e. whether 
it is accreted. This process has b een observed 



in SF XTs such as IGR J 17544-26 19 (iRampy et al 



2009h. as well as in "classical" SGXBs such as VelaX- 
1 dKrevkenbohm et aU 12008 ). At least twice dur- 
ing 02a (at f = 12 and f - 36 ks), the count rate 
from IGR J16207-5129 remained low for several hun- 
dred seconds followed by flares and variable emission. 
These low states are also seen in the hard band (15- 
60keV). Combining the spectra from low-intensity 
states throughout 02, we did not detect a significant 
increase in absorption. Focusing solely on a 2.5-ks 
epoch beginning 1 2 ks after the start of 02, the spectral 
fit yields A^H= (19!^)xl022 

cm , which is compatible 
with the A^H measured outside this epoch. Therefore, 



an occultation by a large wind clump along the line of 
sight to the compact object can not be responsible for 
the dip in luminosity since there is no accompanying 
increase in the A^h- 

The reduction in intensit y could be due to the or- 
bital geometry of the svstem. iPrat et al.l(l2008h showed 
that in die inchned binary IGR J19140H-0951, the A^h 
can be related to the orbital phase reaching a maximum 
at superior conjunction and a minimum at inferior con- 
junction. The hard X-ray flux (20-lOOkeV), which 
is unaffected by photoelectric absorption, was actually 
suppressed at inferior conjunction. This is because the 
shock regions near the neutron star's surface, where 
high-velocity stellar winds encounter the ionized gas 
surrounding the CO's magnetosphere, are partially oc- 
culted by the CO at inferior conjunction. The mini- 
mum flux reported for IGR J 1 9 1 40+095 1 was listed as 
unabsorbed, so this occultation is unrelated to column 
density variations (or the lack thereof). On the other 
hand, if we are truly observing IGR J16207-5129 at 
inferior conjunction, then, as in IGR J19140+0951, we 
should expect the column density at superior conjunc- 
tion to be larger (by an order of magnitude) since the 
X-rays would have to travel through more obscuring 
material. This would imply A^h ~ 10^"* cm"^ at supe- 
rior conjunction. This level of absorption has never 
been observed in IGR J 16207 -5 129 and it is only seen 
in extreme cases such as IGR JI63I8-4848. While 
we can not exclude that this effect is responsible for 
the diminished luminosity of IGR J 16207 -5 129 dur- 
ing 02b, we believe it is unrealistic given the implica- 
tions for the column density at superior conjunction. 

Magnetic and hydrodynamic effects (and possible 
photo-ionization) could be responsibl e for the appar- 
ent drop in luminosity ( BlondinI 1994 ). Movement of 
the CO (and in the case of neutron stars, its ionizing 
magnetic field) through the accretion wind can disrupt 
and even inhibit the accretion of material leading to 
variability and even a temporary attenuation of the X- 
ray source. Such tidal-stream disruptions are accom- 
panied by an increase in the column density of up to 
20-30 times the average density of the wind, which 
we do not observe in IGR J16207-5129. 

The X-ray source may have been partially or fully- 
eclipsed by its supergiant companion star We re- 
mark that a potential ingress phase is noticeable in the 
XIS light curve between MJD 54527.4 and 54527.6 
(Fig E). During the eclipse of such systems, the con- 
tinuum emission is believed to consist of a greater pro- 
portion of scattered X-rays. This should be visible 
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Fig. 6. - — Solutions to a model of an eclipsing SGXB 
(Eq.|6) iRappaport & Joss 1983 ) viewed under differ- 
ent inclination angles (degrees). The CO is assumed 
to be a neutron star of 1.4 Mq, and the companion 
star's mass is set to 20 with a radius of 20 Rq. 
The eclipse is taken to last >30ks (solid curve). The 
dot-dashed lines represent t he limits of orbital radi i 
of clumpy winds in SGXBs (INegueruela et al.luOOSb . 
The shaded region shows the extent of the parameter 
space in which IGR J16207-5129 might be located. 



as an evolving hardness ratio, and it is generally ac- 
companied by an increase in the intensity of the iron 
line relative to th e continuum emission (e.g. VelaX-1: 
Sato et al. and references therein). Count rates 

in the soft and especially the hard bands were low dur- 
ing 02b, plaguing the hardness ratio of this epoch with 
large uncertainties, which made them difficult to com- 
pare with the hardness ratios from 02a. 

In summary, the observed attenuation of the X- 
ray flux can not be explained by an increase in the 
column density since the measured absorption value 
does not change in a significant way between 02a and 
02b, nor between the low and high-intensity states 
of 02. Finally, we mention that whatever mech- 
anism is involved in driving the sporadic emission 
in SFXTs could be responsible for the prolonged 
dip that we see in IGR J16207-5129. The lumi- 
nosity that we measured during the lowest states 
of 02b (2 X 10''"^erg s"') is similar to the quies- 
cent flux reported for the S FXT IGR J 16479 -45 14 
dWalter & Zurita HerasI l2007b . This is an order of 
magnitude greater than the quiescent luminosity ex- 
pec ted from tr ansient X-ray sources (10^""^^ ergs"': 



never been reported for this object implying a high 
duty cycle that is consistent with its classification 
among SGXBs. 

Hence, we are not able to provide any firm conclu- 
sions as to why the observed flux is consistently low 
for such a long time. Nevertheless, the eclipse inter- 
pretation remains a viable explanation and it allows us 
to consider the orbital configuration of the X-ray bi- 
nary. An eclipse duration (At) of at least 30 ks (0.35 d) 
would set the lower limit on the orbital period at 0.7 d. 
We do not see an egress phase in the light curve so 
30 ks should be considered as the lower limit on the 
eclipse duration. Given that a continuous 50-ks XMM- 
Newton observation, which also includes periods of 
low activity (but notably, none are longer than ~2 ks), 
shows no evidence for an eclipse, the lower limit on 
the orbital period should be higher unless the CO is in 
an inclined orbit. 



Following the procedure described in iBozzo et al, 



( l2008h . we can estimate the range of orbital periods 
{F) by using the fact that the separation of the centers 
of mass of the system (a) can be ex pressed as a func- 
tion of the eclipse ha lf-angle 0e via (IRappaport & Joss 
19831: lNagaselll989h : 



(6) 



sin" ; sin" ©„ 



where / is the inclination angle. For ©e, we take ad- 
vantage of a simple trigonometric equivalence: 



20e 

2n 



At 



(7) 



We can then derive the orbital period (related to the 
binary separation a through the generalized form of 
Kepler's Law) for various inclination angles. Figure|6] 
presents the results of this model. Also shown are the 
clumpy winds that are expected in a radius ran ge of 1 .5 



to 2.5 Ri, in SGXBs ( INegueruela et alj|2008h . Based 



e.g. lCampana et al.,,2002) . Such long extinctions have 



on our assumptions, we find that the parameter space 
is bounded by orbital periods between 4 and 9 d with 
inclination angles ; > 50° (shaded region in Fig.|6l). 
The full orbit is within the upper limit of the clumpy- 
wind radius (a < 2.5/?*), unless the orbital period is 
greater than ~10d, which is unusual for SGXBs but 
plausible (e.g. IGR J19140H-0951 has an orbital pe- 
riod of ~14d). For comparison, inclination angles for 
eclipsing HM XBs are gener ally between 60 and 85° 
with ©e ~ 25° ( lNagaselll989l) . 
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5. Summary & Conclusions 

Our Suzaku observation of IGRJ16207-5129 con- 
firms the large intrinsic column density (A^h = 
(16.2;^j j) X 10^~cm"^) that was reported previously, 
while demonstrating the long-term variability of the 
absorbing column. The broadband spectrum covered 
an energy range between 0.5 and 60keV with unprece- 
dented spectral sensitivity above 15 keV enabling us to 
constrain the cutoff energy for the first time. This cut- 
off is situated around 19^^keV which is expected if 
the compact object is a neutron star, although we were 
not able to detect a coherent pulsation period or cy- 
clotron lines. Intriguingly, the observation ends with 
a prolonged (~ 30 ks) state of low activity where the 
luminosity decreases by a factor of 10 on average, and 
up to a factor of 20 between extremes, with little or no 
change in the absorbing column. While our data do 
not permit us to make firm conclusions to the nature 
of this event, one possibility is that the X-ray source is 
being eclipsed by its supergiant companion. 

Confirmation of the putative eclipse in 
IGR J16207-5129 requires further observations of 
the source. Only 4 IGRs are confirmed to be eclips- 
ing binaries: IGR J16418-4532, IGR J16479-4514, 
IGRJ17252-3616, and IGR J18027-2016. The first 
two are classified as SFXTs and the latter two are clas- 
sical SGXBs. The demarcation between SFXTs and 
SGXBs is not as clear as originally believed, and a few 
systems (e.g. IGR J 16479 -45 14, 1GRJ18483-0311, 
and IGR J19294H-1816) might represent intermediate 
states between the SFXT and SGXB classes because 
of their peculiar properties. In this respect, the study 
of IGR J 16207-5 129 will play an integral role in help- 
ing us understand the extent to which the emissivity 
differences between both populations stem from their 
unequal wind and orbital characteristics. 
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